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ABSTRACT 

The TET family of FE(II) and 2-oxoglutarate- 
dependent enzymes (Tet1/2/3) promote DNA 
demethylation by converting 5-methylcytosine 
to 5-hydroxymethylcytosine (5hmC), which they 
further oxidize into 5-formylcytosine and 5- 
carboxylcytosine. Tet1 is robustly expressed in 
mouse embryonic stem cells (mESCs) and has 
been implicated in mESC maintenance. Here we 
demonstrate that, unlike genetic deletion, RNAi- 
mediated depletion of Tet1 in mESCs led to a sig- 
nificant reduction in 5hmC and loss of mESC 
identity. The differentiation phenotype due to Tet1 
depletion positively correlated with the extent of 
5hmC loss. Meta-analyses of genomic data sets 
suggested interaction between Tet1 and leukemia 
inhibitory factor (LIF) signaling. LIF signaling is 
known to promote self-renewal and pluripotency in 
mESCs partly by opposing MAPK/ERK-mediated 
differentiation. Withdrawal of LIF leads to differenti- 
ation of mESCs. We discovered that Tet1 depletion 
impaired LIF-dependent Stat3- mediated gene acti- 
vation by affecting Stat3's ability to bind to its 
target sites on chromatin. Nanog overexpression 
or inhibition of MAPK/ERK signaling, both known 
to maintain mESCs in the absence of LIF, rescued 



Tet1 depletion, further supporting the dependence 
of LIF/Stat3 signaling on Tet1. These data support 
the conclusion that analysis of mESCs in the hours/ 
days immediately following efficient Tet1 depletion 
reveals TetVs normal physiological role in maintain- 
ing the pluripotent state that may be subject to 
homeostatic compensation in genetic models. 

INTRODUCTION 

Embryonic stem cells (ESCs) appear to have an unique 
epigenetic state that maintains the pluripotent genome in 
a stable program of self-renewal, while allowing rapid 
induction of alternate transcriptional programs to 
initiate differentiation (1-8). DNA methylation is one of 
the principle regulators of the epigenetic landscape that 
shapes and refines gene expression programs during 
embryogenesis and stem cell differentiation. While mech- 
anisms of establishment and maintenance of DNA methy- 
lation by DNA methyltransferases are well characterized 
(9,10), it was less clear which enzymatic machinery is 
responsible for DNA demethylation or even which 
pathways lead to DNA demethylation (11,12). Recently, 
it was discovered that proteins of the Tet family of FE(II) 
and 2-oxoglutarate-dependent enzymes (Tet 1/2/3) oxidize 
5-methylcytosine (5mC) to 5-hydroxymethylcytosine 
(5hmC) (13-17), which they further oxidize into 5-formyl 
cytosine and 5-carboxylcytosine, thereby promoting DNA 
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demethylation (18-23). Tet proteins and 5hmC have been 
shown to have roles in cancer and stem cell biology 
(13,16,24,25). Mouse embryonic stem cells (mESCs) defi- 
cient in DNA methyl transferases Dnmtl, Dnmt3a and 
Dnmt3b lack 5mC as well as 5hmC implying that 5hmC 
marks are probably derived from pre-existing 5mC marks 
(26,27). 

In mESCs, Tetl is expressed at high levels, comparable 
to those of the master pluripotency factor Oct4 (16,28). 
Tet2, though well expressed in mESCs, is 5-fold less 
abundant than Tetl (28). Tet 3 is variably expressed in 
many tissues but not in mESCs (16). Though various 
tissues express one or more Tet proteins, the 5hmC modi- 
fication is particularly abundant in mESCs and Purkinje 
neurons (13,14,16,29). In mESCs, Tetl binding and 5hmC 
occupancy are correlated with CpG density, are enriched 
at the promoters and gene bodies of nearly two-thirds of 
all genes, and have been linked to both gene activation 
and repression (27,30-33). Tetl has also been shown to 
be required for the recruitment of transcriptional repres- 
sors Ezh2 (30) and Sin3a (27) at CpG-rich promoters of 
developmental regulators. Though Tetl -mediated regula- 
tion is believed to be due to its catalytic activity, Tetl 
might have other functions in addition to converting 
5mC to 5hmC (27). 

A recent study reported that RNAi-mediated knock- 
down (KD) of Tetl in mESCs resulted in downregulation 
of pluripotency marker Nanog and loss of undifferentiated 
state implicating Tetl in mESC maintenance (16). In 
contrast, subsequent studies found that Tetl KD cells 
were morphologically indistinguishable from control 
mESCs with no changes in Nanog expression, modest re- 
duction in total 5hmC levels, and a minor-to-moderate 
increase in 5mC levels (26-28,31). Discrepancy between 
these studies may be due to differences in KD efficiency, 
off-target effects, or homeostatic compensation masked by 
antibiotic selection. More recently, it was reported that 
Tetl _/_ mice are viable, fertile and grossly normal, and 
that Tetl~ l_ mESCs maintained a normal undifferenti- 
ated mESC morphology with only ~35% reduction in 
total 5hmC levels (34). Since only a modest reduction in 
5hmC levels is observed in this genetic model, it is possible 
that maternal contribution or unaltered Tet2 in Tetl null 
cells may be compensating Tetl's function, thereby 
obscuring the direct impact of Tetl loss. The discrepancies 
observed in the various RNAi-mediated depletion experi- 
ments in addition to the puzzling retention of 5hmC in the 
Tetl null mESCs suggest that unanswered questions 
remain regarding the role of Tetl in the maintenance of 
the pluripotent state of mESCs. 

To understand and clarify the role of Tetl and 5hmC in 
epigenetic and transcriptional regulation of mESCs, we 
used RNAi to acutely deplete Tetl in mESCs and per- 
formed expression and genome-wide 5hmC occupancy 
studies. We find that acute short-term depletion as 
opposed to genetic deletion of Tetl results in a significant 
decrease in 5hmC levels, downregulation of pluripotency- 
associated factors, impairment of LIF-dependent Stat3- 
mediated gene activation, and loss of embryonic stem 
cell identity. 



MATERIALS AND METHODS 

Mouse ES cell culture, RNAi and alkaline 
phosphatase staining 

Oct4GiP mESCs were kindly provided by Dr Austin 
Smith. E14Tg2a mESCs were obtained from Mutant 
Mouse Research Resource Centers, and Jl mESCs were 
obtained from ATCC. The cells were routinely maintained 
on gelatin-coated plates in the ESGRO complete plus 
clonal grade medium (Millipore), and were used at 
passage 20-35 for experiments. For siRNA transfections, 
mESCs were cultured on gelatin-coated plates in 
Ml 5 medium: Dulbecco's Modified Eagle Medium 
(Invitrogen) supplemented with 15% Fetal Bovine Serum 
(FBS) (Invitrogen), 10|iM 2-mercaptoethanol, 0.1 mM 
non-essential amino acids (Invitrogen), lx EmbryoMax 
nucleosides (Millipore), 1000 U of ESGRO (Millipore). 
For cells to be harvested or stained 96 h after transfection, 
20-25 x 10 3 mESCs were transfected with siRNAs at 
100 nM in Ml 5 medium in one well of a 96-well plate. 
About 0.7 |il of Lipofectamine 2000 (Invitrogen) was 
pre-mixed with 10|il of Opti-MEM (Invitrogen) and then 
mixed with 1 x 10 _11 mol siRNAs. Dissociated mESCs 
were plated at 25 x 10 3 per well in 100 |il of Ml 5 medium 
in gelatin-coated plates. siRNA-lipids complexes were 
incubated at room temperature for 15-30min and then 
added to the cells. About 25-50% of the cells in each 
well was re-plated into one well of a 12-well plate the 
next day, and cultured in Ml 5 medium for another 3 
days. For cells to be harvested 48 h after transfection, 
75 x 10 3 mESCs were transfected with siRNAs at 50 nM 
in Ml 5 medium in one well of a 24-well plate, and cells 
were collected 2 days after transfection. RNAi experiments 
were performed using indicated individual siRNAs: Tetl 
siRNA #1 (Invitrogen, MSS284895), Tetl siRNA #2 
(Invitrogen, MSS284897), Tetl siRNA #3 (Dharmacon, 
D-062861-01), Tetl siRNA #4 (Dharmacon, 
D-062861-02), Tetl siRNA #5 (Dharmacon, D-062861- 
03), Tetl siRNA #6 (Dharmacon, D-062861-04) and 
control siRNA duplex targeting firefly luciferase 
(Dharmacon, 5'-CGTACGCGGAATACTTCGA). For 
lineage marker [quantitative reverse-transcriptase PCR 
(qRT-PCR)] and microarray analysis, cells were harvested 
48 and 96 h after transfection for RNA extraction. For 
alkaline phosphatase (AP) staining, cells were fixed, 
permeablized and stained for AP activity with the 
Alkaline Phosphatase Detection Kit (Millipore). 

Nanog overexpression cells 

Mouse Nanog was cloned into pcDNA-EF-HA/Flag 
vector and sequence verified (see Supplementary 
Material). The resulting plasmid, mNanog-pcDNA- 
EF-HA/Flag, was transfected into Oct4GiP cells with 
Lipofectamine 2000 (Invitrogen), and the cells were 
selected with G418 and plated at colonal density. Single 
clones were picked and expanded, and clone #2 was 
chosen for the experiment. The expression of the exogen- 
ous HA/Flag-tagged Nanog was estimated to be at 
the same level of the endogenous Nanog based on 
western blot. 
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Oct4GiP reporter assay 

Oct4GiP cells were transfected in 96-well plates: 0.35 |il of 
Lipofectamine 2000 (Invitrogen) was pre-mixed with 10 \l\ 
of Opti-MEM (Invitrogen) and then mixed with 
5 x 10 _12 mol siRNAs. Dissociated cells were plated 
at lOxlO 3 per well in 100 |il of Ml 5 medium in 
gelatin-coated plates. siRNA-lipid complexes were 
incubated at room temperature for 15-30min and then 
added to the cells. Cells were cultured in Ml 5 medium 
with daily medium change. Four days after transfection, 
cells were lifted and dissociated by trypsinization, pipetted 
into single cell suspension in 10% FBS in Phosphate 
Buffered Saline (PBS), and fluorescence-activated cell 
sorting (FACS) analyzed on LSRII FACS analyzer 
(BD). Percentage of differentiation was determined by 
measuring the percentage of green fluorescent protein 
(GFP)-negative cells. For each experimental condition, 
3-4 independent transfections were carried out, and data 
was plotted as mean ± SEM (standard error of the mean). 

Western blot 

Cells were directly lysed in SDS-PAGE sample buffer. 
Proteins were resolved by SDS-PAGE, transferred to 
nitrocellulose membrane using iBlot (Invitrogen), and 
probed with the indicated antibodies. Loading was 
normalized based on Ran or Tubulin. Primary antibodies 
used: Nanog (Millipore, AB9220), Stat3-Y705-phospho- 
specific (Cell Signaling, 9145), Stat3 (Cell Signaling, 9132), 
Dnmt3a (Santa Cruz, H-295), Dnmt3b (IMEGENEX, 
IMG-184A), Ran (BD, 610340) and Tubulin (Santa 
Cruz, sc-9104). 

Quantitative RT-PCR 

Total RNAs were prepared from cells using the RNeasy 
(Qiagen), and cDNAs were generated using the 
Omniscript RT kit (Qiagen) according to the manufactur- 
er's instructions. Quantitative PCRs were performed using 
the Ssofast Evagreen Supermix (Bio-Rad) on the Bio-Rad 
CFX-384 or CFX-96 Real-Time PCR System. At least 
three biological repeats were carried out for each experi- 
ment. For every biological repeat, triplicate or quadrupli- 
cate PCR reactions were performed for each sample, and 
Actin was used for normalization. Data was plotted as 
mean ± SEM. Please refer to Supplementary Table SI 
for gene specific primers used for RT-PCR analysis. 

5hmC slot blot, and genome-wide mapping and 
analysis of 5hmC 

Mouse E14Tg2a cells were transfected with control 
siRNA or an siRNA targeting Tetl in 6-well plates. 
Cells were lysed 96 h after transfection in Trizol reagent, 
and both total RNA and genomic DNA were prepared 
according to manufacturer's instructions. Total RNA 
was used to determine the Tetl KD efficiency. Genomic 
DNA was quantified with NanoDrop, and 100 or 500 ng 
was loaded on GE Hybond + nylon membrane using the 
GE slot blot device according to manufacturer's instruc- 
tions and blotted for 5hmC. Primary antibody used: 
5hmC (Active Motif, 39769). Relative signal strengths of 



5hmC bands were quantified using the Image-J software 
(h ttp : / /imagej . nih . go v/ij /) . 

For genome-wide mapping of 5hmC, the genomic DNA 
obtained from control KD and Tetl KD cells were 
processed using the QUEST 5hmC detection kit (Zymo 
Research) according to manufacturer's instructions. 
Briefly, the genomic DNA was treated with 5hmC 
glucosyltransferase enzyme, which specifically tags the 
5hmC in DNA with a glucose moiety yielding a 
modified base, glucosyl-5hmC. Next the DNA was 
digested using glucosyl-5hmC sensitive restriction endo- 
nuclease Mspl. Mspl efficiently cleaves the DNA when 
cytosine, 5mC, or 5hmC but not gluey syl-5hmC is 
located at the inner C position within its recognition 
sequence (i.e. CCGG). In other words, the 5hmC sites in 
genomic DNA treated with 5-glycosyltransferase enzyme 
are protected from Mspl digestion. By comparing the 
qPCR or next-generation sequencing results for treated 
versus untreated genomic DNA, 5hmC sites are 
detected (see below for details). Computational analysis 
indicate that digesting mouse genomic DNA with 
Mspl, selecting 40-220 bp fragments, and performing 
36-bp sequencing would cover ~1 million distinct CpG 
dinucleotides (~5% of all CpGs), with roughly half 
located within 'CpG islands' (including sequences from 
90% of all CpG islands) and the rest distributed 
between other relatively CpG-poor sequence features 
(35). Every Mspl cleaved sequence read will thus include 
at least one informative CpG position. Based on this in- 
formation, Mspl-digested DNA from treated and untreat- 
ed samples derived from control and Tetl KD mESCs 
were size-selected (40-220 bp) and sequenced using 
Illumina GAIL 

Sequenced 36-bp reads were aligned to the reference 
genome (mouse NCBI36/mm8 assembly) using Bowtie 
0.12.2 (36). Those reads that mapped to unique genomic 
locations with at most two mismatches were retained. 
Roughly ~ 50-6 5% of all reads were mapped to unique 
genomic locations. To increase the mapping efficiency, we 
trimmed one base from the 3'-end of the unmapped reads, 
aligned the residual 36-bp reads to the reference genome 
and retained those that mapped to unique genomic loca- 
tions. This procedure was repeated until the residual reads 
were of length 18, and, in each iteration, reads that 
mapped to unique genomic locations were retained for 
further analysis. The reads that mapped to unique 
genomic locations but whose corresponding 5'- to 
3 / -sequence in the reference genome neither started with 
CGG nor ended with GGC were considered noise and 
discarded. The rationale for this is because the Mspl 
cleaves the DNA containing its recognition site CCGG 
at the inner C position (C/CGG) (37), the DNA fragments 
that is sequenced must all start with CGG or end with 
GGC (complementary strand). Reads that start or end 
with anything other than CGG or GGC, respectively, 
are likely to be byproducts of DNA degradation or 
sequencing errors. The resulting mapped reads for the 
four samples (control KD treated, control KD untreated, 
Tetl KD treated and Tetl KD untreated) were normalized 
by total reads in each sample so that the number of reads 
mapping to a specific site is directly comparable across all 
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four samples. Sites with less than two reads in all four 
samples were excluded from the analysis. A site was 
defined as a 5hmC site only if the normalized read count 
at this site in the control untreated sample is at least 
1.5-fold higher than that in the control treated sample. 
A 5hmC site was defined to have reduced 5hmC levels in 
Tetl KD mESCs only if 

yvruntreated / yytreated 
A5hmC = iV c°ntrolKD/ jV controlKD > j 5 
A/untreated / A/treated — 
iV TetlKD / iV TetlKD 

where A 5 hmC measures the reduction in 5hmC levels in 
Tetl KD cells and N[ s is the normalized read count in 
sample s with treatment t. Genes were considered to 
have 5hmC occupancy if they have one or more detected 
5hmC sites within 5 kb upstream of the transcription start 
site (TSS) to the transcription end site (TES). Please refer 
to Supplementary Table S2 for PCR primers used for 
quantitative 5hmC analysis. 

Chromatin immunoprecipitation 

E14Tg2a mESCs were crosslinked with 1% formaldehyde 
in PBS for lOmin, and the resulting chromatin was 
sonicated using Diagenode Bioruptor for 15-18 cycles, 
30 s per cycle with a 30 s rest between each cycle, to 
obtain ~200-500bp fragments. Chromatin immunopre- 
cipitation (ChIP) was performed by incubating indicated 
antibodies with Dynabeads protein A magnetic beads 
(Invitrogen, 100-02D) for 4-6 h at 4°C, followed by over- 
night incubation with the chromatin at 4°C. Subsequent 
washing and reverse crosslinking steps were performed as 
described earlier (38). Each ChIP was performed using 
chromatin from ~2x 10 7 cells. Primary antibodies used: 
Stat3 (Santa Cruz, C-20), H3K9me3 (Abeam, ab8898) 
and IgG (Santa Cruz, H-270). Please refer to 
Supplementary Table S3 for ChIP primers used to 
analyze Stat3-target/binding sites. 

ChlP-Seq factor binding data analysis 

Published ChlP-Seq data sets for various transcription fac- 
tors and chromatin remodelers (listed in Supplementary 
Table S4) were downloaded from the Gene Expression 
Omnibus (GEO) (http://www.ncbi.nlm.nih.gov/geo/) or 
the BROAD Institute website. In cases where the data 
was mapped to mouse reference genome assembly mm9, 
the data was realigned to mm8 assembly using the liftOver 
tool (39). Where available for download, binding sites 
(peaks) as defined in the original publication were used. 
Otherwise SISSRs (40) was used for peak calling. Direct 
target genes of a particular factor (for example, Tetl) were 
defined as genes that have one or more detected factor 
binding sites within 5-kb upstream of the TSS to the TES. 

Microarray gene expression data analysis 

Raw CEL files (Affymetrix array type Mouse 430 
2.0) generated from this study were processed using 
relevant R/Bioconductor packages (41). Specifically, 
RMA (42) was used for background correction and data 
normalization. Probe sets were defined using Entrez gene 
based custom chip definition files (CDFs), version 13 



(http://brainarray.mbni.med.umich.edu/Brainarray/Data 
base/CustomCDF/genomic_curated_CDF.asp) (43). This 
resulted in a single expression intensity measure per Entrez 
gene ID. All subsequent analyses were carried out on the 
log 2 scale. To determine differentially expressed genes, a 
family-wise moderated Mest (44) (KD versus control) was 
performed followed by a multiple testing correction pro- 
cedure to control the false discovery rate (FDR) (45). 
Unless otherwise noted, genes were considered differen- 
tially expressed if their FDR < 0.1 and fold change >1.5. 
Since global gene expression changes due to Tetl KDs 
using Tetl siRNAs #1 and #2 were highly similar, we 
treated them as biological replicates in the differential ex- 
pression analysis. Unsupervised hierarchical clustering 
was performed using the CLEAN software package (46). 
Functional enrichment analysis of differentially regulated 
lists of genes was performed using R packages CLEAN 
(46) and GO.db (41). 

Published microarray data sets (listed in Supplementary 
Table S5) were downloaded from GEO (http://www.ncbi 
.nlm.nih.gov/geo/). Where available, raw CEL files were 
re-processed using the above described microarray pro- 
cessing pipeline. Otherwise, the respective Series Matrix 
File was used to obtain preprocessed expression values 
and corresponding probe annotations. All data were 
transformed and analyzed on the log 2 scale. Each probe 
identifier was mapped to an Entrez gene ID using relevant 
Bioconductor packages by continually matching annota- 
tions in the following order until an Entrez gene ID was 
found: Entrez gene ID, RefSeq ID, gene symbol or alias, 
gene name, unigene ID, Ensembl transcript ID. Probes 
mapping to the same Entrez gene ID were averaged 
(after the log-transformation). Fold changes and 
adjusted P-values (KD versus control) using moderated 
^-statistic described above were then computed for each 
data set independently. 

Data availability 

All microarray and 5hmC sequencing data generated for 
this study were deposited in NCBI GEO repository under 
the accession number GSE34267. 

RESULTS 

Acute depletion of Tetl results in loss of embryonic 
stem cell identity 

To elucidate the role of Tetl and 5hmC in epigenetic and 
transcriptional regulation of mESC self-renewal and 
pluripotency, we used RNAi to acutely deplete Tetl in 
Oct4GiP mESCs, which express the GFP under the 
control of the Oct4 promoter (47). Since Oct4 is exclusive- 
ly expressed in ESCs, the GFP expression faithfully cor- 
relates with the ESC identity. Two small-interfering RNAs 
(siRNAs #1 and #2) targeting Tetl were used to ensure 
that the effects were specific. Cells transfected with 
siRNAs were maintained in normal mESC culture condi- 
tions for 4 days, and subjected to FACS to determine the 
fraction of differentiated (GFP-negative) cells. Acute de- 
pletion of Tetl resulted in small but significant increase in 
the percentage of differentiated cells (Figure 1A). 
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Tetl knockdown (KD) efficiency was confirmed using 
qRT-PCR and western blot (Figure IB and 
Supplementary Figure SI A). Examination of Tetl KD 
mESCs for colony morphology and AP staining revealed 
morphological changes and loss of AP staining consistent 
with differentiation (Figure 1C). The outcome remained 
similar even when E14Tg2a or Jl cell lines were used 
(Supplementary Figure SIB). To rule out the possibility 
that the observed changes were due to siRNA off-target 
effects, we used four additional siRNAs (#3-6), some 
of which were used in other studies that observed 
no obvious phenotype (28,31), and confirmed that acute 
depletion of Tetl indeed induces morphological changes 
and loss of AP staining consistent with differentiation 
(Supplementary Figure SIB). 

ESC markers are downregulated and differentiation 
markers are upregulated upon Tetl depletion 

To determine whether Tetl depletion has an effect on gene 
expression programs influencing mESC self-renewal, 
pluripotency and differentiation, we used microarrays to 
profile global gene expression in mESCs transfected with 
control and Tetl siRNAs at 96 h after transfection. 
Pluripotency-associated genes such as Nanog, Esrrb, 
Tell, Tbx3, Klf2, Klf4, Lefty 1, Lefty 2, Tcfcp2ll and 
Prdml4 were downregulated, and differentiation- 
associated genes including ectoderm and neuronal 
markers Fgf5, Pitx2, Nestin, Nefm, CD133 (Proml), 
CD44, Lefl and Zicl and trophectoderm markers Eomes 
and Krt8/18/19 were upregulated in response to Tetl KD 
(Figure ID). Gene ontology analysis of transcripts 
downregulated at 96 h revealed enrichment of many 
terms related to stem cell maintenance and proliferation, 
whereas upregulated transcripts were enriched for terms 
related to development and differentiation (Figure IE). To 
understand the kinetics of downstream effects of Tetl KD, 
we also carried out microarray analysis at 48 h after 
siRNA transfection when 7^7-depleted cells (Figure IB 
and Supplementary Figure SI A) were morphologically in- 
distinguishable from control cells (Supplementary Figure 
SID). Interestingly, our 48 h Tetl KD data clustered 
tightly with those from shRNA-mediated Tetl KD 
studies that did not observe morphological changes even 
past 96 h (Figure ID and Supplementary Figure SIC), and 
away from our 96 h Tetl KD data that accompanies mor- 
phological changes, suggesting a kinetic effect. The lack of 
phenotype in shRNA-mediated depletion studies could be 
due to the (i) differences in reaction rates in shRNA- 
versus siRNA-mediated knockdowns as a result of add- 
itional processing time necessary for the shRNA hairpin 
structure to be cleaved into siRNA or (ii) homeostatic 
compensation of TetTs function masked by antibiotic 
selection. 

We performed qRT-PCR and confirmed gene expres- 
sion changes observed in the microarray analysis. Tetl 
depletion decreased the expression of key pluripotency 
factors: Nanog by nearly 3-fold, and Oct4 and Sox2 by 
small, but reproducible and statistically significant, 
~ 1.25-fold. Other pluripotency-associated genes such as 
Klf4, Esrrb, Tell, Tbx3, Prdml4 and Nodal antagonists 



Leftyl and Lefty 2, which are among the earliest to be 
downregulated upon ESC differentiation (48), were 
downregulated by ~2- to 4-fold (Figure IF). Tetl deple- 
tion also increased the expression of several differentiation 
genes including Pitx2, Nestin, Nefm, Lefl and Zicl that 
are largely associated with Ectoderm lineage (Figure IF). 
Furthermore, these changes were consistent across 
Oct4GiP, E14Tg2a and Jl mESCs (Supplementary 
Figure S2A), and remained consistent even when using 
an siRNA that was used in two previous studies that 
observed no obvious phenotype (28,31) (Supplementary 
Figure S2B). Taken together, we concluded that acute 
RNAi-mediated depletion of Tetl results in the loss of 
mESC identity. 

Tetl depletion results in significant reduction in total 
5hmC levels 

Examination of total 5hmC levels in 7^7-depleted mESCs 
using slot blot revealed a significant reduction in 5hmC 
compared to cells transfected with control siRNA 
(Supplementary Figure S2C). The severity of observed 
morphological changes positively correlated with the 
extent of 5hmC loss (Figure 1G). We suspected that the 
downregulation of Tet2 in Tetl -depleted cells (Figure ID), 
also observed in one other report with findings consistent 
with ours (30), might also be a contributor to reduced 
5hmC levels. Quantitative RT-PCR analysis of Tet2 in 
7ef/-depleted mESCs revealed a reproducible and statis- 
tically significant ~ 1.5-fold reduction in Tet2 levels 
(Supplementary Figure S2D). This along with Tetl occu- 
pancy at Tet2 promoter (Supplementary Figure S2E) sug- 
gested regulation of Tet2 by Tetl. Together, these data 
revealed that the acute reduction in 5hmC primarily due 
to Tetl KD and also possibly due to indirect reduction in 
Tet2 levels is probably the cause for the loss of mESC 
identity. 

Tetl-mediated gene regulation is context dependent 

To investigate the role of Tetl and 5hmC in the mainten- 
ance of mESC identity on the molecular level, we examined 
the genes that were differentially expressed upon Tetl de- 
pletion in our microarray analysis. Using a reasonably 
stringent criterion (FDR < 0.1; fold >1.5), we identified 
530 upregulated and 389 downregulated genes 
(Supplementary Table S7). Examination of published 
Tetl occupancy data (27,30) revealed that Tetl binds to 
72% of upregulated and 63% of downregulated genes with 
no statistically significant preference for one group over the 
other, which is similar to the pattern observed with other 
factors known to play diverse regulatory roles in gene regu- 
lation such as the insulator binding transcription factor 
CTCF (Figure 2 A and Supplementary Figure S3 A). No 
correlation between Tetl binding and gene expression 
changes was observed even when analyzing all genes with 
expression above detectable levels (Supplementary Figure 
S3B). These data indicate that Tetl binding alone is not 
predictive of Tetl's function as an activator or a repressor. 

Tetl has been reported to be required for the recruit- 
ment of transcriptional repressors Ezh2 (30) and Sin3a 
(27) at CpG-rich promoters of developmental regulators. 
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Figure 1. Depletion of Tetl and 5hmC levels results in loss of mouse embryonic stem identity. (A) Oct4GiP mESCs were transfected with indicated 
siRNAs in normal ESC medium and cultured for 96 h. The percentage of differentiated cells was determined by measuring the percentage of 
GFP-negative cells using FACS at 96 h after transfection (**P< 0.001; *P<0.01). Error bars represent SEM of three experiments. (B) Relative 
Tetl mRNA level in control and Tetl KD Oct4GiP mESCs 48 h after transfection. Data are normalized to Actin. Error bars represent SEM of three 
experiments < 0.001). Error bars represent SEM of three experiments. (C) AP staining of Oct4GiP mESCs transfected with control 
siRNA, and Tetl siRNAs #1 and #2. Cells were cultured in normal ESC medium, and AP staining was performed 96 h after transfection. 

(continued) 
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We find that genes upregulated upon Tetl KD are 
enriched for binding of polycomb members Suzl2, Ezh2 
and Ring lb (Figure 2A and Supplementary Figure S3 A), 
whereas downregulated genes are enriched for binding 
of Stat3, Klf4, Oct4, Esrrb and Sin3a among others 
(Figure 2A). Since Tetl binds to nearly two-thirds of all 
genes, its co-occupancy with other transcription factors is 
rather high, which could be misinterpreted for its func- 
tional interaction with these factors. When normalized 
by expected frequencies, Tetl's co-occupancy with 
other factors appears insignificant compared to that 
between factors that are known to function together 
(Supplementary Figure S3C). Like esBAF, a mESC 
specific chromatin remodeling complex (49) that facilitates 
both activation and repression (50,51), we propose that 
Tetl probably plays a similar role as a general facilita- 
tor/recruiter of activating and repressing factors in a 
context dependent manner. 



Tetl depletion results in genome-wide reduction in 
5hmC levels 

To determine the extent of 5hmC loss in Tetl KD cells on 
a genomic scale, we used a modified version of methyla- 
tion sensitive restriction enzyme sequencing (52) and 
generated nucleotide resolution 5hmC maps covering 
^90% of sequences located within CpGs islands in 
control and Tetl KD mESCs (see 'Materials and 
Methods' section for details). We identified 57 895 5hmC 
sites in control mESCs, of which 32 652 had at least 
1.5-fold reduction in 5hmC levels in Tetl KD cells 
(Supplementary Table S8). Quantitative PCR was per- 
formed to validate 5hmC sites and reduction in 5hmC 
levels (Supplementary Figure S3D). We find that genes 
that were upregulated or downregulated upon Tetl deple- 
tion are equally likely to contain 5hmC sites and experi- 
ence reduced 5hmC levels (Figure 2B), which is consistent 
with the established dual roles of Tetl and 5hmC in 
both gene activation and repression (27,30,33). These 
data indicate that an acute depletion of Tetl causes sig- 
nificant global reduction in 5hmC levels that were not 
specific to gene classes, and that the directionality of 
gene expression change (up/down) is probably determined 
by the gene context and local transcription factor occu- 
pancy (Figure 2A). 



Meta-analysis of genomic data suggests interaction 
between Tetl and LIF/Stat3 signaling 

We next examined the fate of Tetl regulated genes upon 
KD or knockout of other pluripotency factors. A 
meta-analysis of microarray data from published studies 
revealed that genes that are differentially expressed upon 
Tetl KD underwent similar expression changes upon 
leukemia inhibitory factor (LIF) withdrawal or after de- 
pletion of key pluripotency-associated factors such as 
Esrrb, Tell, Tbx3 and Kl/2/4/5, known to be associated 
with LIF/Stat3 signaling (8,53-55) (Figure 2C), suggesting 
Tetl's possible interaction with or co-regulation of tran- 
scriptional circuit regulated by LIF/Stat3 signaling 
(Figure 2A). In contrast, Tetl -regulated genes underwent 
largely divergent changes upon deletion/depletion of tran- 
scriptional repressors Tcf3 and Sall4 (Figure 2C). 

LIF-dependent Stat3-mediated gene activation is 
dependent on Tetl 

To examine Tetl's possible mechanistic interaction with 
LIF/Stat3 signaling (Figure 2 A and C), we analyzed 
global gene expression changes upon Tetl depletion 
against those after LIF withdrawal. Tetl -depleted and 
LIF-starved mESCs undergo similar changes in the 
global transcriptional profile (Figure 3 A). The set of 
genes co-activated by Tetl and LIF includes a number of 
factors essential for pluripotency, including Esrrb, Tell 
and Tbx3 (56,57). This suggested that Tetl or Tetl- 
dependent 5hmC levels might be required for LIF- 
dependent Stat3-mediated gene activation. To further 
explore this connection, we decided to study Tetl KD 
mESCs 48 h after transfection, when Tetl is depleted but 
the cells appeared normal compared to control transfected 
cells (Figure IB, Supplementary Figure SI A and D). This 
time point was chosen so that we could clearly distinguish 
cause from effect. First, we asked whether LIF/Stat3 signal 
transduction was impaired in Tetl KD mESCs. No signifi- 
cant changes in mRNA levels for Stat 3 or LIF chimeric 
receptor consisting of gpl30 and LIF receptor (LIFR) 
were observed at 48 h (Figure 3B and C). r<^/-depleted 
mESCs at 48 h after transfection had no significant 
changes in total Stat3 or Nanog levels but a modest reduc- 
tion in phosphorylated/activated Stat3 (Figure 3D). 

Next, we examined the ability of Stat3 to bind to its target 
sites on chromatin in the absence of Tetl. ChIP experiment 



Figure 1. Continued 

See Supplementary Figure SIB for AP staining results for siRNA #1 and four additional siRNAs on Oct4GiP, Jl and E14tg2a mESCs. 
(D) Expression fold changes of selected genes upon Tetl KD in mESCs based on microarray analysis performed 48 and 96 h after transfection. 
Fold changes from data generated for this study [Tetl-KD (96 h) #1, #2 and Tetl-KD (48 h) #1] are presented alongside fold changes observed in 
recently published Tetl KD (>96h) or knockout studies (color-coded on the top-right). Each column corresponds to fold changes obtained from an 
individual array computed in relation to its corresponding control array. Observed/reported morphological changes are symbolically indicated at the 
top, with columns ordered based on unsupervised hierarchical clustering. See Supplementary Figure SIC for a comprehensive heatmap. (E) Gene 
ontology analyses of up- and down-regulated genes in Tetl KD cells compared to control cells. Only selected categories are shown. For complete 
lists, see Supplementary Table S6. (F) Relative mRNA levels of selected mESC pluripotency-associated genes and lineage marker genes in control and 
Tetl KD mESCs at 96 h after transfection. The mRNA levels in control cells are set as one. Data are normalized to Actin. Error bars represent SEM 
of three experiments. (G) Scatter plot showing strong positive correlation among relative Tetl mRNA levels, relative total 5hmC levels, morphogical 
changes, and AP staining in control and Tetl KD E14Tg2a mESCs. Each data point corresponds to a siRNA that was used for transfection. The 
y-axis indicates relative Tetl mRNA levels 96 h after KD, and the x-axis represents quantified intensity of 5hmC signal inferred from slot blot (96 h). 
Bottom panel show the corresponding representative morphological changes and AP staining for each siRNA. Error bars represent SEM of data 
from five replicates (different amounts of DNA were spotted) from two independent experiments. 
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Figure 2. Meta-analysis of genomic data sets suggests interaction between Tetl and LIF/Stat3 signaling. (A) Top: Transcription factor occupancy at 
919 genes differentially expressed upon Tetl KD (389 downregulated and 530 upregulated). Genes are represented along the j-axis and factor 
occupancy is denoted by blue bar. Target gene occupancy is defined as factor occupancy within 5-Kb upstream of the gene's TSS and/or within its 
gene body. Bottom: Log-fold enrichment of factor occupancy at up/do wnregulated genes. Red and green histograms denote enrichment in up- and 
down-regulated genes, respectively. *P<0.01 after Bonferroni adjustment for multiple testing. (B) Top: Presence/absence of 5hmC sites (5hmC) in 
control mESCs, and 5hmC sites with > 1.5-fold reduced hydroxymethylation levels (A5hmC) in Tetl KD mESCs (96 h). Bottom: Log-fold enrich- 
ment at down- versus up-regulated genes, respectively. (C) Expression fold changes of 919 genes differentially expressed upon Tetl KD (96 h; first 
column) presented alongside fold changes observed after KD or knockout (KO) of select other pluripotency factors in published reports. Genes are 
represented along the y-axis. Data sets along the x-axis have been ordered based on unsupervised hierarchical clustering of corresponding gene 
expression fold changes. Blue and yellow rectangles highlight Tetl-hlF and Polycomb clusters, respectively. Cells with no/missing data are colored in 
gray (DKO: double KO). 



using Stat3 antibody was carried out at 48 h after Tetl KD. 
Stat3 binding was affected at many gene targets co-bound 
and co-activated by Tetl and Stat 3 (Figure 3E) indicating 
that Tetl is required to facilitate LIF/Stat3 responsiveness 
in mESCs by enabling Stat3 to bind to its target sites on 
chromatin. Since the expression changes due to loss of Stat3 
binding at Stat3-activated genes is not evident at 48 h but 
only at 96 after Tetl KD (Figure 3E), we concluded that the 
loss of Stat3 binding precedes loss of expression. We do not 
rule out the possibility that a modest reduction in 
phosphorylated Stat3 levels (Figure 3D) might also have 
played a role in impaired Stat3 binding. Regardless, the 
LIF-dependent Stat 3 -mediated gene activation in mESCs 
is dependent on Tetl. 

Nanog overexpression or inhibition of MAPK/ERK 
signaling rescues Tetl KD phenotype 

Since mESCs can be propagated without activation of 
Stat3 by imposing constitutive expression of core 



pluripotency factor Nanog (58), we checked if 
Nanog overexpression (Figure 4A) could rescue the differ- 
entiation phenotype exhibited by Tetl KD cells. As 
expected, exogenous Nanog largely rescued the morpho- 
logical changes and AP activity in Tetl KD cells 
(Figure 4B). Also, since mESCs can be propagated 
without Stat3 activation by suppressing prodifferentiative 
MAPK/ERK signaling (59), we tested if mESCs trans- 
fected with Tetl siRNA in 2i medium (which inhibits 
MAPK/ERK and Gsk-3b signaling) blocks differenti- 
ation. We found that while Tetl depletion slightly 
decreased cell viability in 2i medium (data not shown), 
the expression of pluripotency-associated factors Nanog, 
Tell, Klf4 and Lefty2 were largely rescued (Figure 4C). In 
addition, Tetl KD in 2i medium also suppressed 
upregulation of differentiation genes such as Pitx2, 
CD 133, Nefm and Lefl. Together, these data provide 
evidence to further support the dependence of LIF/Stat3 
signaling on Tetl. 
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Figure 3. Dependence of LIF/Stat3 signaling on Tetl. (A) Two-dimensional gene density heatmap depicting global gene expression changes 96 h 
after Tetl KD against changes observed 48 h after LIF withdrawal (51). The intensity of each square represents the number of genes that fall within 
that square. Axes indicate degree of fold change, from nil (middle of axis) to > 1.5-fold (outermost squares). (B) Gene expression fold changes 
observed in microarray analysis 48 and 96 h after Tetl KD in mESCs, and 48 h after LIF withdrawal (51). (C) Relative Stat 3 mRNA level in control 
and Tetl KD mESCs. The mRNA level in control KD cells is set as one. Data are normalized to Actin. Error bars represent SEM of three 
experiments. (D) Western blot analysis showing protein levels in control and Tetl KD (48 h) mESCs. Tubulin is used as a loading control. (E) Left: 
ChIP assay of select Stat3 target regions using an antibody against Stat3 in control KD and Tetl KD mESCs (48 h). The y-axis represents 
enrichment over input normalized to a negative control region (Yipf2). Refer to Supplementary Figure S4 for data from control ChIP using a 
non-specific antibody against IgG. Error bars represent SEM of three experiments. Right: Gene expression fold changes observed in microarray 
analysis 48 and 96 h after Tetl KD in mESCs, and 48 h after LIF withdrawal (51). 
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Figure 4. Nanog overexpression or suppression of MAPK/ERK signaling rescues Tetl KD phenotype. (A) Western blot analysis showing Nanog 
overexpression with HA-tag. (B) AP staining of E14Tg2a mESCs, with and without Nanog overexpression, transfected with control siRNA and Tetl 
siRNA #1. Cells were cultured in normal ESC medium, and AP staining was performed 96 h after transfection. (C) Relative mRNA levels of selected 
mESC pluripotency genes and differentiation marker genes in control and Tetl KD mESCs in 2i medium. Oct4GiP cells were transfected with 
control siRNA or Tetl siRNA #1 at 50 nM in 24-well plates in 2i-medium (which inhibits MAPK/ERK and Gsk-3b signaling) and cells were 
harvested 96 h after transfection. The mRNA levels in control cells are set as one. Data are normalized to Actin. Error bars represent SEM of three 
experiments. 
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Figure 5. Tetl negatively regulates de novo DNA methyltransferase Dnmt3b. (A) ChIP assay of select Stat3 target regions using an antibody against 
H3K9me3 in control and Tetl KD mESCs (48 h). The j-axis represents enrichment over input normalized to a positive control region (Myodl) for 
H3K9me3. Error bars represent SEM of three experiments. (B) Relative mRNA levels of DNA methyltransferases Dnmtl, Dnmt3a and Dnmt3b in 
control and Tetl KD mESCs 96 h after transfection. The mRNA levels in control cells are set as 1. Data are normalized to Actin. Error bars 
represent SEM of three experiments. (C) Western blot analysis showing protein levels of Nanog, Dnmt3a and Dnmt3b in control and Tetl KD 
mESCs 96 h after transfection. Ran is used as a loading control. (D) Genome browser shot showing a region containing the Dnmt3b gene and results 
from Tetl and Sin3a ChlP-Seq experiments by various groups (GSM numbers denote GEO accession). The red open rectangle highlights Tetl 
occupancy at the promoter region of Dnmt3b, where a CpG island is present (green-filled rectangle). (E) Relative 5hmC levels at Dnmt3b locus in 
control and Tetl KD (96 h) mESCs. Error bars represent SEM of three experiments. 
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Figure 6. Proposed model for 7>/7-mediated epigenetic and transcriptional regulation of mESC self-renewal and pluripotency. Red arrows denote 
regulatory interactions inferred from data generated for this study. Tetl, regulated by Oct4 (28), regulates DNA methylation (5mC) by converting 
5mC to 5hmC (16,28). Tetl regulates LIF/Stat3 signaling by facilitating Stat3 binding by an yet to be determined mechanism, and regulates the 
transcriptional regulatory module comprising Nanog, Esrrb, Tell, Tbx3, Klf2/4, Prdml4 and Lefty 1/2. Tetl's regulation of Tet2 confers tight 
regulation of 5mC to 5hmC conversion. Tetl's negative regulation of de novo DNA methyltransferase Dnmt3b may provide an additional layer 
of Tetl -mediated regulation of 5mC. Silent and active promoters on the chromatin are denoted by broad red and green arrows, respectively. 



Tetl negatively regulates de novo DNA methyltransferase 
Dnmt3b 

Given that genes downregulated in Tetl KD cells undergo 
a minor-to-moderate increase in 5mC levels (26-28,31), 
and that repressive histone mark H3K9me3 has been 
known to recruit DNA methyltransferases and be asso- 
ciated with DNA methylation (35,60), we hypothesized 
that the loss of Stat3 binding and Stat3-mediated activa- 
tion of its targets accompanies an increase in H3K9me3. 
ChIP using an antibody against H3K9me3 revealed that 
this was indeed the case (Figure 5A). It is not clear, 
however, whether Tet proteins alleviate/prevent gene 
silencing at Stat3 targets by preventing deposition of 
H3K9me3, or by maintaining a DNA-hypomethylated 
state. Interestingly, we noticed a 2-fold increase in 
Dnmt3b levels in Tetl KD cells in microarray analysis 
(Figure ID). Quantitative RT-PCR and western blot con- 
firmed the increase in Dnmt3b at the mRNA and protein 
levels (Figure 5B and C). Since Tetl targets Dnmt3b 
(Figure 5D), we examined the 5hmC levels at the 
Dnmt3b locus and found > 3-fold reduction in 5hmC in 
Tetl KD cells (Figure 5E). Based on this evidence, it is 
possible that TetTs regulation of DNA methylation might 
also involve its negative regulation, direct or indirect, of de 
novo DNA methyltransferase Dnmt3b. 



DISCUSSION 

Our studies indicate that acute short-term depletion as 
opposed to long-term genetic deletion of Tetl in mESCs 
results in the loss of the mESC identity. Tetl depletion, in 
our hands, reduced total 5hmC levels by roughly 2- to 
3-fold, Nanog expression by ~3-fold and Tet2 expression 
by ~1. 5-fold. The extent of 5hmC reduction positively 
correlated with morphological changes and extent of loss 
of AP staining (Figure 1) suggesting that it is the loss of 
5hmC and not necessarily Tetl that caused differentiation. 
Our findings are consistent with those by another group 
that used shRNAs to deplete Tetl and observed similar 
changes in morphology, global gene expression, Nanog, 



Tet2 and total 5hmC levels (16,30). Our findings, 
however, are in contrast to Tetl depletion studies by 
others (26-28,31), who observed only ~ 1.5-fold reduction 
in total 5hmC levels, no reduction in Nanog or Tet2 levels 
and no distinguishable morphological changes. Most re- 
markably, our results were reproducible even when using a 
siRNA that was employed in two previous studies that 
observed no obvious phenotype (28,31). Interestingly, 
the molecular changes we observe after Tetl depletion 
are consistent with those from siRNA-mediated KD but 
not shRNA-mediated KD with puromycin selection from 
a previous study (31) (data from Xu et al in Figure ID). 
It should be noted that the only other study that reported 
a Tetl KD phenotype used shRNA-mediated depletion of 
Tetl with no antibiotic selection. Since the siRNAs we 
used do not lead to greater KD efficiency than the 
shRNAs used by studies with no phenotype, it is 
possible that the lack of phenotype observed by others 
could be the result of homeostatic compensation of 
TetTs function masked by antibiotic selection. 

At present, although we have no ready explanation for 
the discrepancy between our results and those from Tetl 
null mice (34), a modest reduction in 5hmC levels in the 
genetic model raises an intriguing possibility of maternal 
contribution or unaltered Tet2 expression compensating 
for Tetfs function after prolonged loss of Tetl, there- 
fore masking the direct effect of Tetl loss that we 
observe hours/days after RNAi-mediated acute Tetl KD. 
Our findings suggest that an acute as opposed to long-term 
chronic reduction in Tetl and thus 5hmC levels is probably 
difficult for mESCs to overcome or compensate, thereby 
resulting in the loss of the pluripotent state. This is 
consistent with Tet2's critical role in regulating 
self-renewal, proliferation, and differentiation in hemato- 
poietic stem cells (61-63), and Tet3's requirement for 
paternal genome conversion of 5mC into 5hmC (64-66). 

The sharp reduction in Tetl and 5hmC impaired 
LIF-dependent Stat3 -mediated gene activation by affecting 
Stat3's ability to bind to its target sites on chromatin, 
indicating dependence of LIF/Stat3 signaling on 
Tetl and/or 5hmC. Suppression of prodifferentiative 
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MAPK/ERK signaling or Nanog overexpression, known 
to maintain mESCs in the absence of LIF or Stat3 activa- 
tion, rescued Tetl KD phenotype further supporting the 
dependence of LIF/Stat3 signaling on Tetl. Collectively, 
our results support a model in which the Tetl -mediated 
regulation of 5hmC levels is critical to the maintenance 
of the pluripotent state in mESCs (Figure 6). Tetl aids 
LIF-mediated Stat3 gene activation by facilitating Stat3 
binding to chromatin by a yet to be determined mechan- 
ism. Tetl regulates key pluripotency-associated factors 
including Nanog, Klf2/4, Tell, Tbx3, Esrrb, Lefty 1/2, 
Prdml4 and Tcfcp211. Regulation of DNA methylation 
by Tetl -mediated conversion of 5mC to 5hmC also 
involves TetFs regulation of Tet2 and Dnmt3b. Our 
findings suggest that studying mESCs in the hours/days 
following depletion of Tetl, before any compensatory al- 
terations could take place, reveals TetFs normal physio- 
logical role in maintaining the pluripotent state, which may 
be subject to homeostatic compensation in genetic models. 
Tetl/2 double knockout may shed light into whether Tet2 
is in fact compensating for Tetl's loss in Tetl null cells. 
Further investigations are needed to better understand the 
mechanistic role of Tetl and 5hmC in facilitating 
pluripotency, and LIF/Stat3 signaling. 
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Supplementary Data are available at NAR Online: 
Supplementary Figures 1-4, Supplementary Tables 1-8, 
Supplementary Methods and Supplementary References 
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